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Abstract 

In the presence of a magnetic impurity the spin-up and down band states are 
modified differently by the impurity. If the multi-electron scalar product (MESP) 
between the occupied spin-up and down states approaches zero then this defines an 
orthogonality catastrophe. In the present paper the MESP is investigated for the FAIR 
(Friedel Artificial Iserted Resonance) solution for a Friedel-Anderson impurity. A 
basis of Wilson states is used. The MESP is numerically determined for the (enforced) 
magnetic, the singlet, and the triplet states as a function of the number N of Wilson 
states. The magnetic and the triplet state show an exponentially decreasing MESP 
as a function of N. Surprisingly it is not the number of states which causes this 
decrease. It is instead the energy separation of the highest occupied state from the 
Fermi energy which determines the reduction of the MESP. In the singlet state the 
ground-state requires a finite MESP to optimize its energy. As a consequence there is 
no orthogonality catastrophe. The MESP approaches a saturation value as function of 
N. 

PACS: 75.20.Hr, 71.23.An, 71.27.+a 
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1 Introduction 



The orthogonality (or infrared) catastrophe was introduced and discussed aheady 40 years 
ago [1], [2], [3j, [4j. An example is a magnetic impurity in a metal host which interacts 
with the conduction electron in the form H' = 2 J (r) s ■ S. The effect of the z-component 
2 J (r) SzSz is the following. Let the spin direction of the impurity point upwards. Then the 
wave function of the conduction electrons is pulled towards or pushed away from the impurity, 
depending on the electron spin. As a consequence the scalar product of corresponding s- 
electron states with opposite spin is slightly less than 1. 

If we denote the resulting (modified) bases for spin up and down as and 

with states in each basis (1 < < A^) and if the lowest A^/2 states are occupied then the 
value of the multi-electron scalar product (MESP) between all occupied s-states with spin 
up and those with spin down is defined by the determinant 



t \J \ /J \J 



Cl,+ |Ci_) \ '^1,4-1 



,t ij \ /J ij 



%/2,+ l'^l,-/ \%/2,+ l%/2- 



The common argument is that the multi-electron scalar product between all occupied 
s-states with spin up and those with spin down approaches zero when the number of oc- 
cupied s-electron states approaches zero. If one compliments this system of impurity with 
spin up plus polarized conduction electrons with the time reversed system where all spin 
directions are reversed then a transition between the two by spin-flip processes of the form 
J (r) [s"'"S'^ + 5^5*^] has vanishing amplitude. Therefore the non-diagonal part of the s ■ S 
interaction tries to prevent the orthogonality catastrophe. This can be well traced in the 
Fair treatment of the Kondo impurity. 

In the following the Friedel- Anderson (FA) impurity will be discussed where this process 
is less obvious. The Hamiltonian for the FA-impurity is given by 



N N 



During the past few years the author has introduced a new numerical approach to the 
Kondo and the FA-impurity problem, the FAIR-method (Friedel Artificially Inserted Reso- 
nance) [S], [S] [Zl- It is based on the fact the n-electron ground state of the Friedel Hamiltonian 
(consisting of an electron band and a d-resonance) can be exactly expressed as the sum of 
two Slater states [8] 

n—1 n—1 

^Fr = Aal n a% + Bd^ [1 (2) 



n—1 n—1 

\Y{a\^^ + Bd^Y{^ 

i=l i=l 



where aj is an artificial Friedel resonance state which determines uniquely the full orthonor- 
mal basis l*^!!- extension of this ground state to the Friedel- Anderson and Kondo 
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impurity problem yields good numerical results. Recently this method was applied to calcu- 
late the Kondo polarization cloud for those impurities [H]. 

Three different solutions of the FA-Hamiltonian will be discussed: the magnetic state, 
the singlet state and the triplet state. For sufficiently large U this Hamiltonian yields a 
magnetic state at temperature only above the Kondo temperature Tk- However, a magnetic 
state can be enforced by the structure of the variational state. This state will be called 
the enforced magnetic state. This avoids the finite temperature treatment. This magnetic 
solution MS has the form 



n— 1 11—1 



MS 



^a!)_i«o+T + + C'^-i^ + ^44 n 4-i n 4+t'^o (3) 



i=l 1=1 



The coefficients A,B,C,D and the compositions of the FAIR states aQ_^ and aj_ are opti- 
mized to minimize the energy expectation value of the FA Hamiltonian. Due to the condition 

(^al^^o \Ho\ aj^$o^ = for i,j>0 and r = +, — the FAIR states determine the other states 

a]^ of the basis jctLl uniquely. 

The singlet state is a symmetric superposition of a magnetic state and its time- (or spin-) 
reversed state, while the triplet state is the asymmetric superposition. (The FAIR states 
and aj_ and the coefficients A, B, C, D are independently optimized for the magnetic, 
singlet, and triplet states). 

For the Fair solution the MESP between the occupied spin up and spin down sub-bands 
is essentially given by 

/ N/2~l N/2-1 \ 

Mf/^) = ( n «l<^oi^oi n s-'^o) (4) 

\ i=0 j=0 I 



2 Numerical Calculation of the Multi-Electron Scalar 
Product 

2.1 The enforced magnetic state 

For most of the numerical calculations Wilson states are used (see appendix A). In the 
calculation the following parameters are used: \Vsd\ = 0.1, U = 1, Ed = —0.5. The magnetic 
solution is optimized for different numbers of Wilson states with = 20, 30, 40, 50, 60. Table 
I shows M{f/^^ of the magnetic solution for the different sizes of the bases. The third 
column gives the scalar product of the two FAIR states, (ao+$oko-'^'o)A/5 5 the fourth column 
the ground-state energy, and the fifth column gives the magnetic moment. ($0 is the vacuum 
state). As one can see the scalar product (ao+$oko-'^'o)A/s ; the ground-state energy (in the 
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enforced magnetic state), and the moment have reached their final values already for N — 30. 
However, the multi-scalar product decreases with increasing N. 



N 










< fl, T\T /o fl, Ar /o ) 
\'^N/2\'^N/2/ J^q 


10 


0.878 


0.823 


-0.607799 


0.514 


.92(.25) 


20 


0.396 


0.501 


-0.627446 


0.687 


.64(.65) 


30 


0.190 


0.4852 


-0.62810 


0.690 


.55(.69) 


40 


0.0917 


0.4845 


-0.62812 


0.690 


.52(.76) 


50 


0.0443 


0.4845 


-0.62812 


0.690 


.50(.77) 


60 


0.0216 


0.484 


-0.62812 


0.690 


.49(.79) 


2*20 


0.394 


0.526 


-0.629323 


0.66 


.64(.61) 


2*30 


0.198 


0.514 


-0.629779 


0.67 


.57(.70) 



Table I: The multi-electron scalar product (MESP) and other parameters for the 
Friedel- Anderson impurity in the enforced magnetic state. The different columns 
give the number of Wilson states, the MESP, the (single electron) scalar product 
(ao+$o|oo-$o)MS between the two FAIR states, the ground-state energy, and the 
magnetic moment. The 6th column is explained in the text. The parameters used 
in the calculation are \Vsd\ — 0.1, U — l,Ed — —0.5. 



In Fig.l the logarithm of the multi-electron scalar product In ^Af^^^^is plotted versus 
the number of Wilson states N. It follows a straight line which corresponds to the relation 

M^^P = 1. 7e-°-°^^*^ = 1.7* 0.93-^ 

Obviously, the multi-scalar product decreases exponentially with increasing N. 




Fig.l: The logarithm of the multi-electron scalar product MESP 

nr=o^ ^I+^o| nr=o^ aJ_$o^ is plotted versus the number of Wilson 
states N with n — N/2 for the magnetic state. 
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In the next step I check whether it is just the number of states N which reduces M}^^ ' . 
For this purpose the N energy cells for = 20 and 30 are sub-divided into two. This is 
achieved by using A = \/2. This doubles the number of Wilson states but adds only one state 
(for positive and negative energy) closer to the Fermi level. The results of this calculation 
are collected at = 2 * 20 and 2 * 30. It turns out that the doubling has essentially only a 
minor effect on m\^!.'^\ This is on a first sight rather surprising since it was believed that 
the increase of the number of states causes the orthogonality catastrophe of the MESP. 

To further confirm this observation I take the energy frame with A^ = 20 and subdivide 
the energy range (—1 : —1/4) into cells with a width of 1/8, replacing two Wilson states 
by six new states. (The same is done for the positive range). This changes M^g from 
0.396 to 0.405. Splitting the same energy range into 14 cells with a width of 1/32 yields the 
MESP M^^g — 0.411. This shows that increasing A^ by subdividing an energy range does 
not contribute to an orthogonality catastrophe (as long as the energy range does not border 
the Fermi level at the energy 0). 

On the other hand, the smallest (absolute) energies have a great impact on the MESP. 
To investigate this question further I take the energies for A^ = 20 and shift the two states 
which are closest to the Fermi level towards the Fermi level. The four energy cells which 
are closest to the Fermi level are = (-2"^ : -2^^^), Cio = (-2"^ : 0), Ui = (0 : 2"^^), 
^12 = (2^^ : 2~^). I replace ±2~^ by ±2~^^. Then the (average) energies of the corresponding 
states are Sg = - j^^ -1. 954 1 x lO^^, ew = -2-^°, en = 2-^^ and £12 = 1. 954 1 x lO'^. 

1 /2 

Of course this reduces the s-d interaction strength Vsd (^) ^ov v = 10, 11 from [2~^/2] = 2~^ 
to [2-^72]^/^ = 2-1°. (The ratio |Kd(z^)|^ /Su remains constant). 

After optimizing the {aj_|_} and {aj_} bases and the coefficients A, B, C, D the resulting 
MESP is reduced to Mj^g = 0.0208. The number of states is still A^ = 20. The shifting 
of the smallest energies from ±2-1° to ±2-^° changes the value of the MESP from 0.396 to 
0.0208. This shows that the value of the MESP is determined by the occupied state closest 
to the Fermi level. The total number of states is only important when it determines the 
energy of this state. 



2.2 The singlet state 

In the next step I calculate the MESP for the singlet ground state. The same parameters 
\Vsd\ = 0.1, U = 1, Ed = —0.5 are used as in table I and Fig.l. The FAIR solution for the 
singlet state is obtained by reversing all spins in ^ms and combining the two states. 

*55 = *M5 (Ti) + *M5 (it) 



5 



n—1 n—1 



+ 



Aal_^al^^ + Bd\al^^ + Cal_^d\ + £>44j 11 "^-i 11 4+i^o (5) 

i=i i=i 

n—1 n—1 



1=1 1=1 



The coefficients A, B, C, D, A', B', C, D' and the compositions of the FAIR states aj_,_ and 
Oq- again optimized to minimize the energy expectation value of the FA Hamiltonian. 
In table II are the corresponding data collected. Again the first four columns give the same 
data as in table I, i.e. the number of Wilson states, the MESP, the (single electron) scalar 
product between the two FAIR states al_^ and al_ and the ground-state energy. The 5th 
column gives the Kondo energy (difference between the relaxed triplet energy and the singlet 
ground-state energy). 

Dependence on the number of states N 



N 




{C'0+\C'0-)sS 


Eo,ss 


AE 


>.999 


10 


0.749 


0.645 


-0.62272 


14.8 X 10-3 




20 


0.742 


0.6448 


-0.637535 


10.1 X 10~3 


9-10 


30 


0.742 


0.6448 


-0.637965 


9.87 X 10-3 


9-21 


40 


0.742 


0.6448 


-0.63798 


9.86 X 10"3 


9-31 


50 


0.742 


0.6448 


-0.63798 


9. 86 X 10-3 


9-41 


60 


0.742 


0.6448 


-0.637973 


9. 85 X 10-3 


9-51 


2*20 


0.751 


0.657 


-0.639684 


10.4 X 10-3 


17-23 


2*30 


0.751 


0.657 


-0.639993 


10.2 X 10-3 


18-43 



Table II: The multi-electron scalar product (MESP) and other parameters for 
the Friedel-Anderson impurity in the singlet state. The different columns give 

the number of Wilson states A^, the MESP, the (single electron) scalar product 
(ao+$oko-*^o)Afs between the two FAIR states, the ground-state energy, and the 
Kondo energy. The 6th column is explained in the text. The parameters used 
in the calculation are \Vsd\ = O.l, U = 1, Ed ~ —0.5. 



For the last column I calculated the scalar product ^al|_ ,j$o|flLj$o/ ^oi all pairs of (i, j) 

which form a A^ x A"-matrix. It turns out that the diagonal elements ^a^|__j$o|al.,i^o^ close to 
the Fermi energy approach the value one. For example if the sixth column shows for A^ = 30 
the value "9 — 21" then the values of the (single particle) scalar products ( aj,. j$o|aLi^o 



lie between 0.999 and 1.000 for 9 < i < 21. Obviously the states a\_,- and ,■ arc almost 
identical in this interval. This is very different for the enforced magnetic state. There, in 
table I the 6th column shows the value of the diagonal scalar product for i = N/2 and the 



larger one of its two neighbors (a]^ j^^^^oW- jv/2±i^o 
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Dependence on the interaction \Vsd\'^ 

The MESP in the singlet state depends on the strength of the s-d interaction. Keeping 
the mimber of Wilson states constant = 40 the MESP is numerically determined and 
collected in table III. 









-£^0,55 


AE 


>.999 


0.10 


0.742 


0.645 


-0.637977 


9. 86 X 10-3 




0.09 


0.708 


0.604 


-0.621392 


8. 13 X 10-3 


10-30 


0.08 


0.665 


0.553 


-0.605078 


6. 14 X 10-3 


10-30 


0.07 


0.607 


0.489 


-0.589200 


4. 09 X 10^3 


11-29 


0.06 


0.527 


0.406 


-0.573975 


2.22 X 10^3 


11-29 


0.05 


0.407 


0.299 


-0.559655 


8.21 X 10-4 


12-28 



Table III: The multi-electron scalar product (MESP) and other parameters for 
the Friedel-Anderson impurity in the singlet state. The columns give the 

s-d interaction IK^I^, the MESP, the scalar product <^aQ^<l>o 1^0-00 ^5 the 
ground-state energy, and Kondo energy AE. The 6th column is explained in 
the text. The parameters used in the calculation are U = 1, E'd = —0.5 and 
the number of Wilson states is A^ = 40. 



In Fig. 2 the logarithm of the Kondo energy is plotted versus the logarithm of the MESP. 
A linear dependence is obtained. The MESP shows a weak dependence on the Kondo energy 
with a power of about 1/4. 




Fig. 2: The log-log plot of the 
Kondo energy versus the MESP 
for the singlet state for different 
s-d interactions. 
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2.3 The triplet state 

If one arranges the spins in each component in equ. ([5]) for the singlet state so that all 
spin-down creation operators are moved to the left and all spin-up ones to the right, then 
the coefficients are pair-wise equal, i.e. A' = A, B' = B,etc. This yields the symmetric or 
singlet state. On the other hand, if one sets the coefficients pair-wise opposite equal, i.e. 
A' = —A, B' = —B,etc then one obtains the asymmetric or triplet state. Of course, one 
has to restart the optimization of ao+, _ and A, B,C, D. For the asymmetric state a 
finite MESP increases the total energy. So if one searches for the relaxed triplet state with 
minimal energy one may expect a strong reduction of the MESP. This is indeed found. In 
table III the data for the triplet state are collected for the same parameters as before. 



N 




(ao+ ao-)T5 


-^0,55 


{(^N/2\(^N/2)rj.f. 


10 


0.926 


0.891 


-0.582822 


.97(.12) 


20 


0.0998 


0.448 


-0.626428 


.16(.85) 


30 


0.0162 


0.481 


-0.628054 


.049(.88) 


40 


2.59x10^3 


0.484 


-0.628117 


.020(.89) 


50 


3.02 X 10-^ 


0.484 


-0.628119 


.082(.96) 


60 


4.17x10^5 


0.484 


-0.628119 


.21(.95) 



Table III: The multi-electron scalar product (MESP) and other parameters for 
the Friedel-Anderson impurity in the triplet state. The columns give the number 

of Wilson states A^, the MESP, ^aj _^$oko -^^o^ s-nd the ground-state energy. 
The last column is explained in the text. The parameters used in the calculation 
are \VJ^ = 0.1, U = 1, = -0.5. 

In the limit of large the MESP approaches zero. That means that the two components 
in equ. ([5]) (top and the bottom line) are completely decoupled. The triplet state consists 
of two orientations of the magnetic state with zero interaction between them. Therefore 
the energies in the magnetic and the triplet state become equal for large A^. A comparison 
between table I and table III does indeed show this agreement. 

The dependence of the logarithm of the MESP on the number A^ of Wilson states is 
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plotted in Fig. 3. 




10 20 30 40 50 60 



Fig.3: The logarithm of the MESP 

(nr="o «l^oi nr="o «l^o) is plotted 

versus the number of Wilson states N 
with n — N/2 for the triplet state 



3 Discussion and Conclusion 



The states aj,_ ^ are constructed from the basis cl by extracting a FAIR state o^^ q. Therefore 



the states a^^_ ^ and cj, are pair-wise quite similar except that there is one state missing in the 
basis 
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Fig.4: The square of the diagonal 



matrix-elements 



aV.,:$o|aL.,;$o 



as well as the ncxt-to-diagonal 



matrix-elements 



are plotted as a function of i for the 
enforced magnetic state. 



As a consequence the two bases j| and quite similar. In Fig.4 the single- 



particle scalar products 



for the enforced magnetic state are plotted as 



a function of i (full circles). In addition 



are plotted as empty up and 



down triangles. The full curve (without symbols) gives the sum of the three contributions. 
One recognizes that an arbitrary state al,_ ^ (for i > 0) can be constructed to 80% out of the 
states j, ^+1 and at On the other hand ^ and at ^ overlap to only 30% for small 
energies (in the center of the horizontal axis) . 
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Fig. 5: The square of the diagonal 



matrix-elements 



as well as the ncxt-to-diagonal 



matrix-elements 



are plotted as a function of i for the 
singlet state. 



This is very different for the singlet state. In Fig. 5 the single-electron scalar products 



$o|al,<I> 



as well as 



$o|al.±i$o 



are plotted as a function of i for the singlet 



state. One recognizes that over a large energy range the states ^ and a[_ ^ are 99% or more 
identical. Only for (absolute) larger energies on the left and right side is the overlap reduced 

to about 70%. 

The reason for this different behavior is rather transparent. The energy expectation value 
of the enforced magnetic state does not depend on the MESP. The s-d transitions happen 
only within the same spin orientation and therefore between the same bases. As an example, 
one has the transition 



n-l 



n-1 



n—1 



n-l 

^i-A n «ii n «m*o <=> 4A n «ii n «It*o 



i=l 



i=l i=l 



Here the matrix element is just 



n-l 



n-l 



n-l 



n-l 



i=l 



1=1 



1=1 



10 



With al_ = the value of (0) is given by 

There are no processes that involve the MESP. 

On the other hand in the singlet state one has transitions from 

tt — l ti—l n—1 n—1 

aj_;aj+^ Yl al_^ aJ+^$o <=> 11 "J-T 11 "I+i^o 

i=l 1=1 1=1 i=l 

Such a transition is proportional to the square of the MESP. To be able to harvest energy 

4. 



from these processes the states aL and aj_ are, for small energies, aligned parallel to each 



other. 



for the 



In Fig. 6 the corresponding single-particle scalar products 

triplet state are plotted as a function of i. With exception of the few (2 - 4) points in 
the center the results is very close to the result of the enforced magnetic state. The two 
(magnetic) components of the triplet state are essentially decoupled because the coupling is 
proportional to [M*^^/^] which is of the order of 10^^. Therefore the triplet state is essen- 
tially equal to the sum of the enforced magnet state plus its time- (spin-) reversed partner 
(for sufficiently large N). 
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Fig. 6: The square of the diagonal 



matrix-elements 



«'+,.'^o|a-,i<^'o 



as well as the ncxt-to diagonal 



matrix-elements 



are plotted clS 9; function of i for the 
triplet state. 



The results of this paper are two-fold. One has to distinguish whether (the expectation 
value of) the energy in the ground state depends on the MESP. If Eq is independent of the 
MESP and one uses the Wilson states as the original basis then the MESP M^l"^ decreases 
exponentially with increasing N . This is due to the fact that the energy of the states 
closest to the Fermi level decreases exponentially as well. Just by moving these states closer 
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decreases M^/^. On the other hand if the total number of states in a given energy interval 
is increased without reducing the energy of the states closest to the Fermi energy then the 
MESP jg bg^^giy affected. 

If the energy in the ground state depends on the MESP (as for the singlet state) then 
this results in a freeze of the MESP. This achieved by forcing the new basis states aj_,_ and 
aj _ to be parallel within a certain energy range of the Fermi level. 



4 Appendix 

A Wilson states 

Wilson [To] in his Kondo paper considered an s-band ranging from —1 to 1 with a constant 
density of states. In the next step Wilson replaced the energy continuum of s-states by a 
discrete set of cells. First the negative energy band is subdivided on a logarithmic scale. The 
discrete energy values are -1, -1/2, -1/22,-2"'^, .. - 2-(^/2-i)^ g. These discrete = -2"^ 
points are used to define a sequence of energy cells: the cell (ti, (for i/<N/2) includes all 
states within {^^-i '■ ^u) = (—1/2'^"-'^ : —1/2'^). A new (Wilson) state cj, is a superposition 
of all states within an energy cell {C,u~i '■ C,u) and has an (averaged) energy + ^u) /'^ = 

^— — . This yields a spectrum 8^. — |, — |, — ••5 — ^jl72) "^f/i- This spectrum is 

extended symmetrically to positive energies (for v > N/2). The essential advantage of the 
Wilson basis is that it has an arbitrarily fine energy spacing at the Fermi energy. For a given 
the two smallest energy cells extend from ±2"*^^/^^^^ to 0, and the (absolute) smallest 
energy levels are ±2~^/^. 
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